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Abstract. The pyrochlore lattice Heisenberg antiferromagnet has a massive 
classical ground state degeneracy. We summarize three approximation schemes, 
valid for large spin length S, to capture the (partial) lifting of this degeneracy 
when zero-point quantum fluctuations are taken into account; all three are related 
to analytic loop expansions. The first is harmonic order spin waves; at this order, 
there remains an infinite manifold of degenerate coUinear ground states, related 
by a gauge-like symmetry. The second is anharmonic (quartic order) spin waves, 
using a self-consistent approximation; the harmonic-order degeneracy is split, but 
(within numerical precision) some degeneracy may remain, with entropy still of 
order L in a system of sites. The third is a large- A' approximation, a standard 
and convenient approach for frustrated antifcrromagnets; however, the large-A'^ 
result contradicts the harmonic order at 0(5") hence must be wrong (for large S). 

PACS numbers: 75.25.-fz,75.10. Jm,75.30.Ds,75.50.Ee 



1. Introduction 

The defining property of a "higlily frustrated" magnet is massive classical ground state 
degeneracies as in quantum Hall systems or Fermi liquids, the high density of (zero 
or) low-energy excitations facilitates a rich variety of correlated states. [T] In three 
dimensions, the pyrochlore antiferromagnet. realized in A2B2O7 oxides or in B sites 
of AB2O4 spinels [2], is considered the most frustrated case [3l [J. We ask what is 
its ground state for quantum Heisenberg spins with large S, till now an unresolved 
question 

In experimental pyrochlore systems, this degeneracy is most often broken by 
secondary interactions (e.g. dipolar [7], Dzyaloshinski-Moriya, or second- neighbor 
exchange) or by magnetoelastic couplings [5l [H [9]. Nevertheless, the pure model 
demands study as the basis for perturbed models, and perhaps to guide the search for 
systems with exceptional degeneracies: Heisenberg models can be cleanly realized by 
cold gases in optical traps [10] . 

The pyrochlore lattice consists of the bond midpoints of a diamond lattice, so the 
spins form corner sharing tetrahedra, each of which is centered by a diamond site. We 
take Ns to be the number of spins (pyrochlore sites) , and L to be the linear dimension. 
We have Heisenberg spins, with 5* ^ 1, and nearest-neighbor couplings Jij = J, so 




(1.1) 
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where Lq, = X^iea ^ tetrahedron spin. (Here a, Ukc other Greek indices, always 
runs over diamond sites, and "i G a" means i is one of the fom' sites in tetrahedron 
a.) The classical ground states are the (very many) states satisfying Lq = for all 
tetrahedra. 

The obvious way to break the degeneracy is the correction energy E' from 
perturbation about some tractable limit, such as: (i) Holstein-Primakoff expansion 
(1/5), as in Sec.[2]and[31 below; (ii) Large- TV expansion, as in SeclU or (iii) expansion 
about the Ising limit of XXZ model [iT]. But which degenerate states to expand 
around? Commonly, one just computes and compares E' for two or three special 
states that have exceptional symmetry or a small magnetic cell. 

Instead, our approach is to express E' as an effective Hamiltonian Ti^^ [12], 
for a generic classical ground state, often via crude approximations that have no 
controlled small parameter, yet result in an elegant form. For any H'^^ , we seek (i) 
its (approximate) analytic form (ii) its energy scale, (iii) which spin pattern gives 
the minimum i?harm, and (iv) how large is the remaining degeneracy. The effective 
Hamiltonian has value beyond the possibility (as here) that it leads us to unexpected 
ground states. First, we can model the T > behavior using a Boltzmann ensemble 
exp{—PH^^) [12]. Second, starting from H'^^ , more complete models may be built 
by the addition of anisotropics, quantum tunneling [T3], or dilution [T^]. Apart from 
analytics, we also pursued the brute-force approach of fitting to a database of 
numerically evaluated energies; minimizing the resulting Ti°^ may well lead us to a 
new ground state not represented in the database. 

Our analytic approach was devised anew for each model; still, a common thread is 
to manipulate the Hamiltonian till the Ising labels of the discrete (collinear) states (see 
below) appear as coefficients in the Hamiltonian, and expand, even though there is no 
small parameter. It is no accident that the effective Hamiltonians are always written in 
terms of loops [5l[TTl[T4] in the lattice, or that the degeneracy-breaking terms have such 
small coefficients. Indeed, all collinear states would be exactly symmetry-equivalent 
if our spins were on the bond- midpoints of a coordination-4 Bethe lattice [T5]fin 
place of the diamond lattice). (This same Bethe lattice will also provide an excellent 
approximation for resumming a subset of longer paths for our loop-expansions.) 

In the rest of this paper, we summarize three calculations [HI [TT] [181 [HI HO] 
for the T = ordered state of the large-S" nearest-neighbor quantum antiferromagnet 
on the pyrochlore lattice. In each real-space expansion produces an effective 

Hamiltonian in terms of products of spins around loops. Sees. [2] and [3] are based 
(respectively) on the harmonic and quartic order terms in the spin-wave expansion. 
In effect, we have a hierarchy of effective Hamiltonians, each of which selects a small 
subset from the previous ground state ensemble, yet still leaves a nontrivial degeneracy 
(entropy of 0{L)). Sec. [His based on large-A^ mean-field theory, an alternative way 
to see anharmonic effects, where the additional limit is taken of a large length for the 
Sp{N) "spin"; the large- loop expansion is different, but gives similar in form 
to the anharmonic case. Along the way (Sees. 12.31 and 13. 3p we comment on related 
models, such as the kagome or checkerboard antiferromagnets, as well as field-induced 
magnetization plateaus. Finally, a conclusion (Sec. [5]) speculates on the prospects to 
address spin-disordered ground states. 
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2. Harmonic effective Hamiltonian 

For sufficiently large S, an ordered state is expectecU as spin fluctuations become sclf- 
consistently small [HI [22] ■ The first way the classical degeneracy may be broken is 
the total zero-point energy of the harmonic spin-wave modes: 

i^harm({ni}) = ^ ^fi^m, (2.1) 
m 

where w„i are the frequencies of all spinwave modes fluctuating around a particular 
classical ground state = Siii with unit vectors {Ai}. (Strictly speaking, the 
constant term in (|2.2ap . below, should also be counted with -Eharm-) This is implicitly 
a function of the local classical directions {fii}, and can be considered an effective 
Hamiltonian that breaks the classical degeneracy. In any exchange-coupled system 
-E'harm IS cxpcctcd to bc a local minimum in collinear [231 124] states, such that all 
spins are aligned along the same axis (call it z), thus hi = rjiz: we assume this from 
now on. From = 0, every tetrahedron a has two up and two down spins. 

2.1. Holstein-Primakoff (l/S) expansion and spinwave modes 

Eq. (j2.ip is the expectation of just one term in the Holstein-Primakoff expansion of 
the Hamiltonian with l/S* as the small parameter: Ti ~ £'ciass — JSNs + Tiharm + 
Tiquart + 0{S~^), whcrc E'ciass is thc classical (mean- field) energy, and 

Hharm = J ^ J ^ CT, . cr, ; (2.2a) 

Wquart = ^^^I'^^i'^] " \^ ^ ' '^J ('^^ + '^j) 

(ij) 

(2.26) 

Here J= J(1+1/2S'); henceforth, we fix J = 1. We choose to expand in spin deviation 
operators cr, = {af,a^), defined so that Oj = (r/icrf +ia\)/^/2S is the standard boson 
operator that lowers the component of spin parallel to n^. The harmonic term 
(|2.2a|) only appears to be independent of the {rji\, which label distinct classical ground 
states; the dependence is hidden in the commutation relations, [uf , crj] = iSrjiSij. The 
anharmonic terms (j2.26| will be the basis for Sec. [3l here the brackets include all four 
combinations of x or y with i or j. 

For any classical ground state, half of the modes are "generic zero modes" [17] 
and have uJm = 0. The other half are "ordinary" modes. Finally, "divergent" modes 
are special ones with divergent fluctuations; these occur where the generic-zero and 
ordinary branches become linearly dependent. It can be shown in real space that 
a divergent mode's support can be bounded to an irregular slab normal to a (100) 
coordinate axis [T7]. Hence the divergent modes have an 0{L) degeneracy, and in 
Fourier space are restricted to lines in (100) directions. The elastic neutron structure 
factor should have sharp features along divergence lines. As we shall see in Sec. [3l 
divergent modes dominate the anharmonic corrections to the energy. 

I This does not contradict the evidence for spin-disordered (spin Uquid or valence bond crystal) 
states at 5 = 1/2 |21l |4]. or in the classical case [3]. 

§ Footnote 13 of Ref. 1161 noted, in the spirit of 1241 . that J^'^m same for any classical 

configuration. But in collinear states, '^ut'^ attains a maximum, which makes it plausible that 
^harm o<: X] l"^! a, minimum. This was confirmed, numerically, for the pyrochlore model in |17| . 
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2.2. Trace expansion and loop effective Hamiltonian 

The {u!m} in (|2.ip are the same as the eigenfrequencies of the (hnearized) classical 
dynamics, [31 \W\ which reduces to 

5La = -S'J^/iQ/jz X (JL^. (2.3) 

P 

This defines an important matrix fi with elements ^at3 = '7i(a^), where i{a, (3) is the 
pyrochlore site that links neighboring diamond lattice sites a and f3; fiap = if a — /3 
or the diamond sites are not neighbors . (In Ref. |17| the same matrix is derived more 
rigorously from the quadratic form in (j2.2Qp .) Thus, via the trick of using tetrahedron 
spins, the dynamical matrix is the classical Ising configuration {rji}. If we can only 
massage the formulation so it appears as a perturbation, an expansion will generate 
the desired effective Hamiltonian. 

The eigenvalues of fx^ are (HuJm/S)^, so the harmonic energy (|2.ip is 

i?ha™({'7j) = J^TrQ[/x2]V2^ , (2.4) 

We can formally Taylor-expand the square root in (|2.4p about a constant matrix Al 
in powers of (/x^/4— ^i)". (Naively, yl = 1 since the diagonal of /x^/4 is the identity; 
actually, larger A is needed to account for additional contributions cx 1 from higher 
powers of fj,.) After collecting powers of fi, we have 

1/2 oo 

= 5^C2feA-(^-i/2)Tr/x2fc(2.5) 
where the coefficients {c2k} have closed expressions. Now, Tr(/x2fc) 

is a sum over all 

of the diagonal terms of /x^*^, i.e. a sum over products of ^t^^ along all of the closed 
paths - on the diamond lattice - with 2k steps. These paths may retrace themselves, 
which gives trivial factors rjf = 1; but steps that go once around a loop contribute 
a configuration-dependent factor ±1 equal to the product of Ising spins around that 
loop. To assure convergence of the sum in (|2.5|) . A> 1.4 is needed. 
Thus, we can re-sum (|2.5|) to obtain an effective Hamiltonian 

H^^^harm = ^^0^. + ^6$6 + A's^S + • • • , (2.6) 

where $2; is the sum over all products Y[ Vi taken around loops (without acute angles) 
of 21 spins in the pyrochlore lattice. 

Most of the retraced path terms are in 1-to-l correspondence with paths on 
the coordination-4 Bethe lattice. This gave a quite accurate approximation for the 
constant term Eq , as well as for the contributions from higher powers of fi that resum 
to give each K21 coefficient in (j2.6p . Then expanding (|2.5p up to the I — 30 term |T7] 
(i.e. loops of length <60), and extrapolating to n = 00, wc obtained the coefficients 
in (EH): Eo = -0.56407VS' , Ke = 0.01365, A'g = -0.0033S'. To test (EH), we 
numerically computed the zero-point energy (j2.ip for many coUinear ground states. 
As confirmed by Fig. [TJa), Ti.'^^hmm represents the energy well. 

2.3. Gauge-like symmetry, ground state degeneracy, and discussion 

The exact harmonic energy admits "gauge-like" transformations, relating one Ising 
configuration to another: rj^^ ^-^ = TaT/}r]^a,f3) where Tq = ±1 arbitrarily on every 
diamond site. In matrix notation, /x' = Tfj,T~^, where r = diag({Ta}). Then /i' is 



S'Tr 



Al 
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Figure 1. Numerical results for degeneracy-breaking energy versus effective 
Hamiltonian written in terms of loop; each point represents a different 
configuration of Ising spins (a) Harmonic spin waves (Sec. [2]l (b). Quartic 

spin waves (SecfSjl; only 7r-flux states (harmonic ground states) are included, (c) 
Large- Af /large-S approximation (Sec. |4j. 



similar to /x and has the same eigenvalue spectrum, so £'harm({'7i}) = ^'harm ({'?;}) ■ 
These are not literally gauge transformations, since the classical ground-state condition 
must independently be satisfied: namely, X^iea = in every tetrahedron. Since 
-E^harm is gaugclike invariant, its value can only depend on gaugelike invariant 
combinations of {jyi}, i.e. loop products, which explains why (j2.6p has exactly the 
form of a Z2 (Ising) lattice gauge action. 

It follows that the harmonic ground states are degenerate; as ()2.6|) implies and the 
numerical calculation confirmed, they are all the (coUinear) configurations in which 
the loop product is Y\-qi = —1 around every hexagon. We call these it -flux states in 
the language of Ref. [6]. The divergent modes of Sec. 12. II provide a trick to construct 
and count gauge-like transformations and hence the ground state degeneracy. Namely, 
any gauge-like transformation can be factorized into two, involving the even and odd 
diamond lattice sites; these transformations in turn correspond one-to-one with a 
basis of divergent modes. In this fashion, an upper bound [T7] on the ground state 
entropy was obtained, of order LlnL. On the other hand, a lower bound of order L 
is easily obtained by explicitly constructing a subfamily of 7r-fiux states by stacking 
independent layers of thickness a/4 in (say) the [001] cubic direction. Each layer is a 
set of chains running in the [110] or [110] direction, with spins alternating both along 
and transverse to the chains, so there is a twofold choice for each layer [TB]. 

Two important loose ends of our harmonic calculation are (i) it was not proven, 
but only checked numerically [5l [17], that coUinear states are local minima of the 
harmonic zero-point energy (j2.ip as a function of classical orientations; nor was it 
proven that they are the only stationary points, (ii). We do not yet understand the 
full set of harmonic (tt-Aux) ground states for the pyrochlore: only a special subset 
are given by the layer stacking construction [T7] . 

The (harmonic-order) loop expansion is easily adapted to similar Heisenberg 
antiferromagnets that support coUinear ground states, such as the checkerboard 
lattice [6]. More interesting are the kagome or pyrochlore antiferromagnet at a large 
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field "magnetization plateau" [TT1[T7|: in that case, the signs get reversed in (|2.6p . (In 
the pyrochlore case, with tetrahedra, H^^^harm now favors a positive loop product 
$6 = +1.) In each case, the ground state entropy again comes out 0{L). 

Hassan and Moessner [^Hj got corresponding results for the W^harm of kagome 
antiferromagnets (including noncoUinear states) in a (variable) field, uncovering 
further subtleties of the degeneracies. Also, Bergman et al [11] extended the derivation 
of (|2.6p to their easy-axis limit. 



3. Anharmonic spin- wave theory 

At harmonic order, the pyrochlore antiferromagnet has wriggled loose from our efforts 
to pin it down to a unique ground state. Evidently, we must try again using the 
anharmonic terms. As in the kagome case, brute-force perturbation theory ~ taking 
the expectation (Tiquart) in the ground state of Tiharm ~ fails, since the harmonic 
fluctuations are divergent. Instead, we must construct a reasonable ground state by 
using the anharmonic terms self-consistcntly. 



3.1. Self- consistent decoupling 



The quartic term Tiquart can be decoupled in a standard fashion: in each quartic 
term, simply pair the operators and replace one of the pairs by its expectation, in every 
possible way. The result is a "mean- field" Hamiltonian - quadratic like Tiharm, but now 
all eigenfrequencies are nonzero (except the Goldstone mode) and all divergent modes 
have been regularized. One can interpret its ground state wavefunction variationally 
as being the best harmonic-oscillator state for the actual Hamiltonian. The effective 
nearest-neighbor interactions are modified as 



J ij ^ Jij ^" ^Jij^ ^Jij 



1 



1 



(3.1) 



Here, 



Gij = 



is the correlation function of fluctuations, which we 



can evaluate numerically. The anharmonic energy depends on a completely different 
set of modes than the harmonic energy did. In light of (|3.1[) . iJquart is dominated by 
the divergent (at harmonic order) modes introduced in Sec. 12.11 those are zero modes, 
which do not contribute iSharm at all (recall ()2.ip ). In principle, then, our recipe is to 
guess a regularized Hamiltonian, compute its correlations {Gy }, and insert these in 
p.ip to get a new Hamiltanian; then, iterate until this converges. 



3.2. Mean field Hamiltonian and self- consistency 

We need to understand the {Gy} due to divergent modes. These modes simultaneously 
enjoy all properties of both ordinary and generic zero modes (since these branches are 
becoming linearly dependent). We use the fact that any ordinary mode [T7], satisfies 

„(«)(,) ^1^^ ^ (3.2) 

where {ui'"''} is the eigenvector of {/^a/s} having the same eigenvalue Wm, and where 
the sum runs over both diamond sites linked through pyrochlore site i. 

Note that as 5 oo, Tiquart ^ Wharm, SO we approach a pure harmonic 
Hamiltonian. In this limit the modes are almost gaugelike-invariant (the regularization 



Effective Hamiltonians for large-S pyrochlore antijerromagnet 



7 



breaks the invariance). We may assume the lower-order term i?harm has been 
minimized, i.e. a 7r-flux state. Any such state is specially uniform in that all hexagons 
are the same, modulo the gauge-like symmetry, and hence all bonds and sites are 
equivalent, Thus, it turns out, the fluctuations of the diamond-site modes in (j3.2p have 
correlations with a simple form parametrized by constants gi,g2, and 173: (w^) = go 
(the same on every site); (uaU^) = riaprjp^g2 for second neighbor (a, 7) on the diamond 
lattice, having (3 as their common neighbor (here g2 < 0); and (by bipartiteness) 
{uaUp) ~ for nearest neighbor (a,/3). Inserting into (|3.2p . we find the correlations 
are Ga = ga and Gy = ^[rnVjdo + 92]- Substituting this into finally gives 

^Ji] = " 2^ (fo + V^Vj\92\) = -SJ* - ^inVr (3-3) 

The constant is absorbed in a re-renormalization to J* = J(l — SJ*); the key (small) 
parameter is e, which breaks the gauge-like invariance and cuts off the divergences. 

Thus, the mean- field Hamiltonian is well approximated by the simple form p.3p . 
which in effect says "strengthen the satisfied bonds relative to the unsatisfied bonds." 
That is also the simplest possible form of a variational Hamiltonian that is consistent 
with the local spin symmetries. In practice, we simply assumed (j3.3[) and computed 
(72(e), so the problem reduces to one self-consistency condition, e/8 ~ \g2{e)\/2S'^ . It 
turns out Gij ^ g2 ~ 5'lne, hence e ^ lnS'/5 and finally iJquart ^ (InS*)^. (The log 
divergence is a consequence of the strongly anisotropic momentum dependence of the 
modes near the divergence lines in reciprocal space.) 

3.3. Effective Hamiltonian, numerical results, and discussion 

We calculated the quartic energy numerically for various periodic states, grouped in 
families within which the states are gaugelike equivalent. They had unit cells ranging 
from 4 to 32 sites, and five gauge families were represented, in particular the tt- 
flux states (ground states of W^harm)- When the result was fitted to an effective 
Hamiltonian, all but a few percent of the anharmonic energy is actually accounted 
by gauge-invariant terms, of the same form as p.6p . The gauge-dependent energy 
differences between states are much smaller in the 7r-flux state, and larger in a gauge 
family where the gSMge- invariant loops are most inhomogeneous. We searched for the 
optimum among (we believe) all possible 7r-flux states in the several unit cells we used 
(with Ns up to 192 sites). 

We performed a numerical fit to an effective Hamiltonian of the form 

H-^^^quart = C^{S)V& + ^8(^)7^8 + CiQ (5)7^10 , (3.4) 

where Vi is equal to the number of loops of length I composed solely of satisfied bonds 
(i.e. with alternating spins). Here 6*6(100) « -0.0621, C8(100)/iV^ « -0.0223. (These 
energies were fitted to (In 5*)^ dependence [18], as implied by the analytics; but our 
range of S values is too small to distinguish from some other power of InS*.) The 
scatter plot in Fig. [Hb) shows that the fit (|3.4p captures the leading order dependence 
on the Ising configuration that splits the harmonic-order degeneracy. 

To explain (j3.4p analytically, note the quartic energy is proportional to the 
energy of (|3.3p . evaluated as if it were harmonic. That can be handled by a small 
generalization of the loop expansion of Sec. O The leading state-dependent term 
turns out to be Vq (with Cg oc e^), confirming analytically a form we had originally 
conjectured empirically. 
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The highest and lowest energy states of p.4p have, respectively, the smallest and 
largest numbers of hexagons with spin pattern tititi EO]- The maximum fraction 
(1/3) of such alternating hexagons is found in a set of states constructed by layering 
two-dimensional slabs. Within our numerical accuracy, these are degenerate for any 
S. These stacked states have the same number of alternating loops of all lengths up to 
16 (we checked), indeed (we believe) up to 26. We conjecture a tiny splitting of these 
states at that high order, maybe even smaller than the similar case of our large-A^ 
calculation (see Sec. 14.31 below.) 

We also calculated the anharmonic effective Hamiltonian for the (planar) 
checkerboard lattice, a tractable test-bed for pyrochlore calculations [3l[6l[26]. But, 
inescapably, two bonds of every "tetrahedron" , (appearing as diagonals of a square) 
have no symmetry reason to be degenerate with the other four bonds. So, in the 
anharmonic calculation, the diagonal bonds renormalize to be weaker than the rest, 
and a unique ordered state is trivially obtained. 

Bergman et al |llj developed a quite different derivation of effective Hamiltonians 
- nicely complementary to ours - by expanding around the Ising limit. Their 7i°^'s 
have a form quite like (|3.4p - i.e., the terms count the number of loops with different 
Ising configurations - and identifying the ground states is comparably difficult. It 
would have been valuable if we had generalized our anharmonic calculation to the 
magnetization-plateau case (see end of Seed]), to compare with the results of Ref. [TT]. 
(The complication of this generalization is that (j3.3[) will get a term with 77,; -I- r]j , 
necessitating a second variational parameter in addition to e.) 

One naively, but wrongly, expected a similar story for the Heisenberg quantum 
antiferromagnet on the pyrochlore lattice as on the (previously studied) kagome 
lattice mi m]. There, spin wave fluctuations selected coplanar (not coUinear) 
configurations as local minima; since all bond angles are 120°, all coplanar states 
were degenerate to harmonic order unlike our result in Sec. 12.31 Due to the non- 
collinearity, the counterpart of (|2.ip had a term Hcubic, third order in {c^/^}, and 
O(^cubic) contributed the same order as O(Hquart) [221 [27]. A consequence was 
distant-neighbor terms in the effective Hamiltonian [57] (which selected the unique 
"•\/3 X -^3" state). A second consequence of non-collinearity was that all generic 
zero modes - an entire branch ~ were divergent, and hence in the kagome case, the 
anharmonic energies (and squared spin fluctuations) both scaled as 0{S^^^), much 
larger than the (In 5)^ of the pyrochlore case. 

4. Large-iV approach to large-S* limit 

Besides the spin-wave expansion, there is another systematic approach to go beyond 
basic mean field theory: Schwinger bosons. Each (generalized) spin has Sp(A^) 
symmetry and is written as a bilinear in boson operators {biam}, where cr =t, i and 
TO runs over TV flavors; the representation is labeled by k which generalizes 25*. The 
physical case is SU(2) = Sp(l), but the N ^ 00 hmit can be solved exactly and is 
often successful as a mean-field theory or the starting point of a 1/N expansion; [28] 
this is popular as an analytic approach to S* = 1/2, in small k limit, since exotic 
disordered ground states can be represented as well as ordered ones. [28] In our work 
[19] . (with Prashant Sharma as the major collaborator, who initiated us into this 
approach) we instead pursued the large- approach for large k. This gives an easier 
recipe for ground state selection than the spin-wave approach, since in large-A the 
degeneracies are usually broken at the lowest order [28] . 
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But which saddle point to expand around? In the pyrochlore antiferromagnet, 
there are exponentially many, corresponding to the same coUinear states as in the 
spin-wave expansion and labeled by the same Ising variables {r)i\. Prior studies just 
investigated high symmetry states, or every state in a small finite system [551 [HI • We 
pursue instead the effective Hamiltonian approach. 

Large-N mean field theory 
Exchange interactions are quadratic in "valence bond" operators Qij = m ~ 

b\l^rJ']^,m^ and in the boson number operator = Y.a.rrMa,,J'i'y,m.- (Here m < N 
is a flavor index in the large generalization, and a j.) For the physical SU{2) 
spins, we have • Sj N^N^ — QljQij- A decoupling quite generally gives 

^5p(JV) = IY. (^IQyf + Q^Aj + H.c) (^^ " (4-1) 

with the classical numbers Qij = {Qij)/N. The Lagrange multipliers Xi, which (it 

turns out) have the same value A = 4k on every site, enforce the physical constraint 

that the boson number is exactly k (the generalized spin length) at every site. We 

want the first nontrivial term in a 1/k (semiclassical) expansion. 

The desired ordered state is a condensation of bosons, (feicr,m) = VNSi^mXia- 

The mean-field ground state energy, obtained via a Bogoliubov diagonalization, is 

^Sp{N) ~ ^fp{N) + ^Sp(N)- the first is the same in every classical ground state. The 

quantum term is the bosons' zero-point energy: 

N , 

EspiN) - ({QyDylTr VA^I - QtQ -iV,A]. (4.2) 

In a coUinear classical state, Qij = K{r]i — '7j)/2, i.e., ±k for every satisfied bond but 
zero for unsatisfied bonds. 

4-2. Loop expansion and effective Hamiltonian 

We have manipulated -E5p(Ar) into the form of a trace of a matrix square root, as in 
(|2.5p for the harmonic spin wave energy - but this matrix Q connects pyrochlore sites 
i, whereas fi in Sec. [2] connected diamond-lattice sites. A Taylor expansion of (|4.2p 
gives the desired effective Hamiltonian, 

^*<».^-f f:|5i^TV(Q.Qr (43) 

rn — 1 

Evidently Tr[(Q^Q/K^)'"] is just the number of closed paths of length 2m on the 
network of satisfied bonds; this network is bipartite, so every nonzero element of 
QtQ is K^. Those paths that eventually retrace every step can be put in 1-to-l 
correspondence with paths on the Bethe lattice (more precisely, a "Husimi cactus" 
graph |19|). They contribute only a constant factor independent of {77^} as do 
paths decorated by additional loops that lie within one tetrahedron. The effective 
Hamiltonian is a real-space expansion in loops made of valence bonds: 

n'^'^spiN) = ^ [CaNs + CeVe + CsA + • • •) (4.4) 

where V21 is the number of non-trivial loops of length 21 with alternating spins, on 
(now) the pyrochlore lattice. The coefficients {C21} were given as a highly convergent 
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infinite sum. hence could can be evaluated to any accuracy: we got Cg = — 3.482 x 10"'^, 
Cg = —3.44 X 10""', and C21+2IC21 ~ 1/10, so short loops dominate. 

Our large- iV loop expansion can be extended to all non-collinear classical ground 
states [ini HO], with the form of (|4.4p but with generalized V2i- It can be applied 
the kagome and checkerboard lattices, [19j giving the usual answers for their ground 
states [281126]. 

^.3. Numerical results and discussion 

We calculated the self-consistent energy for many different collinear classical ground 
states, obtained by a random flipping algorithm described in Ref. [T^. Eq. (|4.4[1 is an 
excellent fit of the state-dependent energy even with just the 21 = & and 2Z = 8 terms, 
as shown in Fig. [Ijc). An independent numerical fit agreed to within 1% for Cg and 
10% for Cg [H]. 

We performed classical Monte Carlo simulations of the Ising model with (|4.4p as 
its Hamiltonian to systematically search for the ground state, using large orthorhombic 
unit cells with 128 to 3456 sites. The optimum was found for a family of nearly 
degenerate states built as a stack of layers, so the entropy of this family is 0{L). 
(Each layer has thickness 3a/4 and there arc four choices per layer, but this is not 
the family found in Sec. O) These states are a subset of those with the maximum 
value Ve — Ns/3 - i.e., one-third of all Ng hexagons are tititi ~ a-nd Vs = 23Ns/6. 
However, it turns out a tiny energy difference ~ 10~^ per spin, corresponding to the 
21 = 16 term in (|4.4[) . splits these states and selects a unique one. 

Let us check our results against the spin-wave approach of Sec. [2l The harmonic 
term of the 1/5 expansion must dominate at sufficiently large S, so the physical 
(SU(2)) semiclassical ground state musthe a ground state of that term, namely a "tt- 
flux" state. Yet the ground states of (j4.4p are not 7r-flux states, and therefore cannot 
possibly be the true ground state: the 1/N expansion has let us down. Nevertheless, 
if Espi^N) values are compared within the "gauge" family of 7r-flux states, the ordering 
of these energies is similar to the quartic spin- wave result (Sec. [3]). 

5. Conclusion 

The trick of writing the zero-point energy as the trace of a matrix, (Eqs. (j2.4p 
and (j4.2p ) - and, for the spinwave expansions, transposing to the diamond lattice 
(Eqs. (j2.3p ) and p.2p ) - enabled an (uncontrolled) expansion giving the effective 
Hamiltonian in terms of Ising spins as a sum of over loops. In each case, there was a 
degenerate or nearly degenerate family of states with entropy of 0{L). The practical 
conclusion is clear, at least : beyond harmonic order, energy differences are ridiculously 
small and would not be observed in experiments. 

Is it, then, possible to realize a disordered superposition of these states, once we 
add to our effective Hamiltonian the "off-diagonal" terms, representing the amplitudes 
for tunneling between collinear states? (Compare [13] for the kagome case, and [TT] 
for the pyrochlore.) Unfortunately, the 0{L) entropy of ground states implies that 
transition from one to another requires flipping 0{L^) spins, so the tunnel amplitude 
is exponentially small as L 00. One also noticed that collinear selection (mentioned 
before Sec. 12. ip provides a different route than "spin ice" to realize an effective Ising 
model in a pyrochlore system; when these collinear states are allowed tunnelings 
(i.e. ring exchanges), won't we realize the "?7(1) spin liquid" of [29]? To stabilize 
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a quantum superposition, the tunnel amplitude should be larger than the energy 
splittings among collinear states, but smaller than the energy favoring coUinearity - yet 
in the pyrochlore, both energy scales are comparable (of harmonic order, i.e. relative 
order 1/S'). The kagomc lattice - or in d = 3, the garnet lattice of corner-sharing 
triangles - is far more promising for disordered spin states, since its harmonic-order 
ground states have extensive entropy. 
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